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ABSTRACT 

Aims. We explore modifications to the current scenario for the slow neutron capture process in asymptotic giant branch (AGB) stars 
to account for the Pb deficiency observed in post-AGB stars of low metallicity ([Fe/H] ^ -1.2) and low initial mass (^1-1.5 Mq) in 
the Large and Small Magellanic Clouds. 

Methods. We calculated the stellar evolution and nucleosynthesis for a 1.3 Mq star with [Fe/H] = -1.3 and tested different amounts 
and distributions of protons leading to the production of the main neutron source within the *^C-pocket and proton ingestion scenarios. 
Results. No ^-process models can fully reproduce the abundance patterns observed in the post-AGB stars. When the Pb production is 
lowered the abundances of the elements between Eu and Pb, such as Er, Yb, W, and Hf, are also lowered to below those observed. 
Conclusions. Neutron-capture processes with neutron densities intermediate between the s and the rapid neutron-capture processes 
may provide a solution to this problem and be a common occurrence in low-mass, low-metallicity AGB stars. 
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1. Introduction 


During the past decade significant information has been gathered 
on the chemical compositions of post-asymptotic giant branch 
(AGB) stars in the Milky Way. This has led to the discovery 
of a class of post-AGB stars that have C/O > 1 and display 
extreme enrichments in the abundances of the elements heav- 
ier than Fe produced by slow neutron captures (the process , 

I Van Winckel & Revniersll200f)t iRevniers & Van Winckell[2003t 
IRevniers et al.ll2004l) . Since AGB stars can become C-rich and 
have been confirmed both theoretically and obser vationally as 
the m ain stellar site for the s process (see, e.g., IBusso et all 
12001 . it is natural to interpret these post-AGB observations 
as the signature of the nucleosynthesis and mixing events that 
occurred during the preceeding AGB phase. These events are 
currently identified as: (i) the mixing of protons into the ra¬ 
diative He-rich intershell leading to the formation of a thin re¬ 
gion rich in the main neutron source (the “pocket”), 
(ii) proton-ingestion episodes (PIEs) inside the convective ther¬ 
mal pulses (TPs), and (iii) the third dredge-up (TDU), which 
carries C and i-process elements from the He-rich intershell to 
the convective envelope and to the stellar surface. Since the de¬ 
tails of all these processes are very uncertain (see discussion 


in, e.g .. [Busso et alJll999l : lHerwidl200^ICamr)bell & Lattanziol 
l2008t) . observations of post-AGB stars provide strong observa¬ 
tional constraints. Recent observations of the chemical com¬ 
position of four low-metallicity ([Fe/H] from -1.15 to -1.34), 
i-process-rich, C-rich post-AGB stars in the Large (J050632, 
J052043, and J053250) and Small (J004441) Magellanic Clouds 
(EMC and SMC, respe ctively) have provided a challenge to 
AGB 5-process models (IDe Smedt et al.l 1201 2[ Ivan Aarle et al.l 
I 20 TI iDe Smedt et al.11^411 ‘L Since we know the distance of 
these stars, from the observed luminosity it is possible to de¬ 
termine that their initial stellar mass was in the range 1-1.5 
Mq. Stellar AGB models in this range of mass and metallicity 
can produce the high observed abundances of the s-process el¬ 
ements, such as Zr and La (1 < [Zr/Fe] < 2 and 1 < [La/Fe] < 
3), together with [Pb/La] ^ 1, if a deep TDU is assumed after 
a last TP. Instead, negative [Pb^ a] values are observed as up¬ 
per limits (iDe Smedt et al.llToidl) . Here we test different possible 
modifications of the current AGB s-process scenario to explain 


' We do not discuss the composition of the m ildly j-process- 
enhanced J053253 reported bv ivan Aarle et alj ( 120131) because this star 
has been recently classified as a young stellar object candidate, indi- 
cating that its abun dances more likely reflect the initial composition 
dKarnath et al.ll201^ . 
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the neutron-capture abundance pattern observed in the MC post- 
AGB stars. 


2. The stellar models 


We have computed the structural evolution of a star of ini¬ 
tial mass 1.3 Mq and [Fe/H] = -1.3 using the version 
of the Mt Stromlo/Monash evolutionary code updated by 
ICampbell & Lattanziol ( |2008|) . Low tempera ture opacities have 
been f urther updated to those calculated by iLederer & Aringei 
(2009). Mass loss was included using the formula of iReimers 
( 1975 ) d uring the RGB pha s e (wit h rj - 0.4) and the for¬ 
mula of IVassiliadis & WoodI (1 19931) during the AGB phase. 
Instantaneous mixing was used in convective zones and the con¬ 
vective bou ndaries were defined usi ng a search for “convective 
neutrality” (iFrost & Lattanzioll 19961) . No overshoot was applied 
beyond this, which resulted in a small number of TDU episodes; 
5 in total - from the 8* TP to the 12* TP - out of a total of 15 
TPs, with the amount of mass carried to the envelope, Mtdu^ of 
the order of a few 10“^ Mq. The last three TPs (13, 14, and 15) 
occurred after the star had left the AGB track (Fig.[TJ. During the 
15* TP the post-AGB star became a born-again AGB star. In our 
model this TP did not experience the TDU nor the ingestion of 
the entire H-rich envelope of mass ~ 0.002 Mq. The observed 
post-AGB stars did not experience such ingestion either since 
they are not H deficient. 

We compute d detailed nucle osynthesis models using a post¬ 
processing code (ICannonlll993h with a network of 320 nuclear 
species from neutrons to ^^°Po. Nuclear reaction rates were 
from the JINA database as of May 2012. The init ial abundances 
in the nucleosynthesis models were taken from lAsplund et alJ 
(l2009l) scaled to the required [Fe/H], except for O, which was 
a-enhanced such that [0/Fe]= -1-0.4 both in the evolutionary and 
in the nucleosynthetic calculations. During the post-processing 
we artificially introduced a given abundance Yp of protons at 
the deepest extent of each TDU episode, which led to the for¬ 
mation of the pocket. Table [1] shows the main features of 
the '^C-pocket models where Yp is the proton abundance, Mpmz 
the mass extent of the p-rich region, and M*' the total inte¬ 
grated mass of protons. We run all of the models twice; with 
protons inserted after each of the five TDU episodes, and with 
protons inserted only after the last TDU episode (cotTespond- 
ing to the 12* TP), except for pocket_case2 which was run with 
protons inserted only after the last TDU episode. As done by 
iDe Smedt et al.l (1201 2l) . we assumed a further last TDU episode 
to artificially occur in all our models after the 13* TP (when the 
mass of the envelope is 0.033 Mq). For this TDU episode we 
set Mtdu to the values given in Table [1] such that the observed 
[La/Fe] is matched for each of the four stars. The TDU is very 
uncertain as it is affected by both physical and num erical choices 
(see, e.g.. lFrost & Lattanzioll 1 9961 ; iMowlavilll 9991) . however, the 
amount of TDU does not affect the abundance patterns discussed 
below but only alters the absolute abundances. 

The proton abundance distribution used in pocket_casel 
(thin line of left panel in Fig.|2]l is the sa me or very simila r to the 
choic e made in the models presented by iDe Smedt et al.l (120121 
12014 and, as expected, produces the same results. According 
to our stellar evolution model the always burns completely 
before the onset of the following TP. In pocket_case2 we experi¬ 
mented with artificially keeping the temperature below 70 MK in 
the pocket, which is not enough to activate the '^C(a',n)'®0 reac¬ 
tion. A much lower value of the rate o f this reaction is excl uded 
since it is known within a factor of ~A (iBisterzo et al.ll20f5h . but 
it is possible that a TDU episode leading to the formation of the 



Fig. 1. HR diagram showing the evolution of our 1.3 Mq, 
[Fe/H] = -1.3 model from the start of the main sequence to the 
white dwarf cooling track. The locations of the four post-AGB 
stars are represented by the red dots with error bars. The loca¬ 
tion of the third last thermal pulse (TP 13) is also indicated. (The 
color version of this figure is available in the online journal.) 


pocket occutTed in the earlier pulses, when the core mass 
and the temperature are lower. In fact, the core mass at which 
the TDU begins is un certain and may be lower than in our model 
(iKamath et alJl20l'^ . For example, if the first TDU episode oc¬ 
curred after the 2"*^ instead of the 8* TP the temperature in the 
region of the pocket would not have exceeed 70 MK before 
the development of the next TP. In this case the neutrons are re¬ 
leased inside the following TP as the is ingested. Next, we 
set the proton abundance Yp to be constant with values equal to 
0.70 X 10“"^ (pocket_case3) or 1.05 x lO^”* (pocket_case4, thick 
line of left panel of Fig. |2l and pocket_case5). Physically, this 
may correspond to efficient mixing inside the pocket dur¬ 
ing the interpulse period induced by the difference in the angular 
momentum between the core and the envelope in a rotating AGB 
star. This mi xing carry the neutro n poison ''^N into the ^^C-rich 
layers fe.g.. IPiersanti et ani2013l) . lowering the s-process effi¬ 
ciency. However, the actual effect is uncertain because magnetic 
fields can slow down the core by coupling it to the envelope, 
and inhibit mixing. In the radiative '^C-pocket pocket_case4 the 
neutron flux lasts for roughly 30,000 yr, the total time-integrated 
neutron flux (neutron exposure) is ^0.5 mbarn ', and the neutron 
density reaches a maximum of 5.2 x 10® cm“^. In the case of the 
convective pocket (pocket_case2) the neutron flux lasts for a 
much shorter time, of the order of 10 yr, and the neutron density 
reaches a much higher maximum of 2.5 x lO" cm^^ for a similar 
neutron exposure. 

We also run some preliminary tests to simulate a proton- 
ingestion episode (PIE) during the post-processing by inserting 
an artificial proton abundance distribution in the intershell before 
the convective zone of TP 13 reaches its maximum extension. At 
this time the H shell is completely extinguished. These protons 
are ingested as the convective region extends outward. As in the 
'^C-pocket models, we artificially simulated a single last TDU 
event mixing the material from the He intershell into the enve¬ 
lope after TP 13 had subsided, setting the amount of TDU mass to 
the value required to match the observed [La/Fe]. Table|2]shows 
the main features of the PIE models, where Y™“ is the maximum 
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Table 1. Properties of the '^C-pocket models. Mtdu is given for the models that include only one pocket. 



Yp 

burning” 

Mpmz(Mq) 

M“(Mo) 


artificial Mtdu( 10 3 M©) 







J004441 

J050632 

J052043 

J053250 

pocket_casel 

standard*’ 

radiative 

1.8 X 10“^ 

1.35 X lO-'* 

45.3 

1.05 

3.13 

4.16 

pocket_case2 

standard!’ 

convective 

1.0 X 10“^ 

7.50 X 10-3 

3.77 

0.18 

0.51 

0.67 

pocket_case3 

0.70 X lO-'* 

radiative 

1.8 X 10-3 

1.26 X 10-3 

87 

1.33 

4.03 

5.40 

pocket_case4 

1.05 X lO-'* 

radiative 

1.8 X 10-3 

1.89 X 10-3 

6.04 

0.27 

0.78 

1.02 

pocket case5 

1.05 X lO-'* 

radiative 

3.6 X 10-3 

3.78 X 10-3 

2.76 

0.14 

0.39 

0.50 


Notes. “Refers to the usual case where all the burns in radiative conditions during the interpulse period or the case where we forced the '^C 
to burn in convective conditions while it is ingested in the following TP. ^’Refers to a proton abundance Yp that decreases exponentially with the 
mass depth from the envelope value to 10“'* at a mass depth Mpmz below the base of the envelope (Fig.O. 



Fig. 2. Proton abundance distribution artificially inserted during the post-processing in our models to simulate the formation of a 
pocket after the 12* TP (left panel) and a PIE event during TP13 (right panel). Left panel: The thin line represents the standard 
exponentially decreasing proton distribution (pocket_case 1) and the thick line the case where the proton abundance Yp is set to 
constant equal 1.05 x 10 (pocket_case4). The dotted line represents the depth reached by the TDU. Right panel: The thin and 
thick lines represent two different proton distributions with a similar total amount of protons ingested; PIE_casel (thin line) and 
PIE_case2 (thick line). The dotted line represents the maximum extention of the intershell convective region associated with TP13. 
The steep proton profile left over by the H-burning shell is also shown in the plot just above 0.638 Mq. 


Table 2. Properties of the PIE models. 



‘Y'max 

P 

M“x.(Mo) 

M“(Mo) 


Mtdu(10 3 M©) 


PIE.casel 

0.035 

8.0 X lO-"* 

4.48 X 10-*^ 

J004441 

0.68 

J050632 

0.03 

J052043 

0.10 

J053250 

0.13 

PIE_case2 

0.007 

3.5 X 10-3 

5.25 X 10-® 

1.08 

0.06 

0.16 

0.20 

PIE case3 

0.007 

4.0 X 10-3 

6.0 X 10-® 

0.41 

0.02 

0.06 

0.08 


proton abundance, Mext the mass extent of the p-rich region, M*' 
the total integrated mass of protons, and Mtdu the TDU mass re¬ 
quired to match each star. The Mtdu required in the PIE models 
are typically lower than in the '^C-pocket models. 

The proton abundance distribution is set to exponentially de¬ 
crease with mass from a value Y™™ = 0.035 (PIE_casel, thin 


line of right panel of Eig. |2]l or 0.007 (PIE_case2, thick line of 
right panel of Eig.|2]and PIE_case3) to a value of 10 over Mext. 
The protons are ingested in the convective region and captured 
by '^C, which results in the production of and subsequently 
of free neutrons for the s process via '^C(a;,n)'®0. We ran a se¬ 
ries of models varying the total mass of protons ingested M*‘ 
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and found best solutions for Mp°^ = 4.5 - 6 x 10 ® Mq 
M ext = 0.8 - 4 X 10“^ Mq, depending on Y™“. This mass extent 
is relatively large: a more realistic profile to simulate overshoot 
of the convective intershell region into the proton tail left over by 
the H-burning shell would have a higher Y”™ (e.g, ~ 0.2) and a 

much smaller Mext (~ 10“^ Mq, see Fig.|2]l. However, to set and 
follow such a profile appropriately requires very high temporal 
and spatial resolution, leading to extremely long running times 
with our computational tools. We leave this task to future work. 

In comparison to the scenario of the pocket, the neutron 
flux lasts for a very short time, of the order of 10 ' yr, five orders 
of magnitudes shorter than in the radiative pocket, while the 
neutron density reaches a maximum of 5 x 10'^ cm~^ (e.g., in 
the PIE_case2 model), seven orders of magnitude higher than in 
the radiative pocket. Even at this value of the neutron den¬ 
sity, the neutron-capture path does not move further away from 
the valley of j3 stability by more than two or three atomic mass 
numbers, for each given element, and still proceeds sequentially 
throughout each value of the proton number. 

Einally, because the PIE was included in the post-processing 
calculations only, we cannot evaluate its feedback on the stellar 
structure and evolution. PIEs can lead, for example, to a splitting 
of the convective region, as energy released by H burning creates 
a temperature inversion inside the convection zone. The amount 
of protons ingested in our model is small (up to 6 x 10 ® Mq) 
but this has been found in some PI E models to be enough to start 
the splitting of the c onvective zone dMiller Bertolami et al.l2005 
iHerwig et al.ll20li1) . Howev er, this effect is mode l dependent, 
e.g., the 3D PIE models of IStancliffe et al.l (1201 ih do not find 
the split to occur e ven for larger amounts of protons, while the 
3D PIE models of I Woodward et al] (1201 5h do. The method we 
have used to simulate a PIE is clearly artificial and preliminary, 
but any ID model not considering hydrodynamics is necessarily 
inaccurate, and also 3D PIE models are not in agreement with 
each other. Nevertheless, the usefulness of our approach lies in 
the opportunity of identifying major qualitative differences be¬ 
tween the PIE and the '^C-pocket scenarios. 

3. Results 

If protons are introduced after each of the five TDU episodes 
computed by the stellar evolution code it is not possible to find 
a solution for any of the four stars since the predicted [Pb/Ee] 
for the models that match [Zr/Ee] is always above the observed 
upper limit (by -1-0.2 dex for J050632 and by -1-0.6 for the other 
three stars). When we include only one pocket, instead, the 
production of Pb is lower and it is possible to find a case that 
simultaneously matches [Zr/Ee] and [Pb/Ee] (Eig. [3]l. The rea¬ 
son for this difference is that with more pockets the material 
is recurrently exposed to several neutron fluxes, which pushes 
the abundance pattern towards Pb. In pocket_case2 Mpmz con¬ 
trols the total amount of ingested and consequently the neu¬ 
tron exposure. The best fit to the observations was found using 
MpMz = 10“^ Mq. a main difference in the results between the 
convective and the radiative models is that the [Rb/Ee] ratio (at 
Z-31) is similar to [Zr/Ee] in the convective case while it is 1 dex 
lower than [Zr/Ee] in the radiative case. This is due to the higher 
neutron density, which favours the activation of the branching 
point at ^^Rb producing th e neutron magic ^^Rb (Ivan Raai et al.l 
l2Q12l:[K^akas et ^l2Q12h . The PIE model results are shown in 
Fig. El Because the general i-process elemental abundance pat¬ 
tern is a one-to-one function of the neutron exposure, it is pos¬ 
sible to obtain very similar results in the PIE models to the '^C- 


pocket models by setting the amount of ingested protons to ob¬ 
tain a similar amount of free neutrons. As in the case of convec¬ 
tive pocket, [Rb/Ee] is close to [Zr/Ee] in the PIE cases due 
to the high neutron density. Observations of this element would 
be a powerful tool to distinguish between the different scenarios, 
however, it is difficult to estimate its abundance in these stars be¬ 
cause it presents only one very strong line, which is saturated. 

The '^C-pocket models typically predict higher [C/Ee] ra¬ 
tios than observed, however, this result depends on Mpmz- For 
example, in pocket_case5 we doubled Mpmz and matched the 
observed [C/Ee], within the error bars, for the three EMC stars. 
To match the [C/Ee] in J004441, instead, requires a tenfold in¬ 
crease of the mass of the ^^C pocket, which would cover the 
whole intershell. Extra-mixing phenomena at the base of the en¬ 
velope during the red giant and AGB phases could remove C to 
produce N. Using the upper limit of the observed C abundance, 
the results for pocket_case5, and assuming that the excess C is 
converted into N via extra mixing, we expect [N/Ee] in J004441 
to be greater than 2.5 dex. This is two or ders of magnitude large r 
than the value of 0.61 + 0.5 derived by iDe Smedt et akl (l2014l) . 
The PIE models, instead, typically predict lower [C/Ee] ratios 
than observed, due to the lower Mtdu- The O abundances are 
generally underestimated by all the models, which indicates that 
a higher O abundance must be present in the intershell and could 
be achieved if the timescale of He burning during the TPs was 
longer or if some mixing occur red between th e base of the con¬ 
vective zone and the C-O core (lHerwig||2000l) . 

The four stars show positive [Y/Zr] while all the models pre¬ 
dict negative values. We do not pursue this issue further because 
there may be a systematic error in the observations since the Y 
abundance is often based on a limited number of spectral lines. 
Another serious problem is that the models that produce [Pb/Ee] 
below the observational upper limits together with Y and Zr rea¬ 
sonably close to the observations underproduce the abundances 
of the elements between Eu and Pb, such as Eu, Dy, Er, Lu, Yb, 
Hf, and W. The number of lines used to determine the abun¬ 
dances of these elements are limited and often blended, however, 
we cannot ascribe the mismatch to observational issues since all 
the four stars clearly show this problem for several of these ele¬ 
ments. 

4. Discussion and conclusions 

It is not possible for the s process to produce high abundances of 
the elements from Eu to W together with a Pb deficiency. This 
is because once the bottleneck at the nuclei with magic num¬ 
ber of neutrons 82 ('^^Ba and ^^®La) is bypassed the neutron- 
capture flux proceeds through each isotope to accord¬ 

ing to its neutron-capture cross section. In other words, the 
abundance pattern between the magic neutron numbers at Ba 
and Pb is almost completely determined by the neutron-capture 
cross se ctions of the isot opes involved. These are relatively well 
known (iBao et alJl2000l) and that their values cannot be drasti¬ 
cally changed is also demostrated by the fact that the i-process 
has no major problems (within 10-20%) in reproducing the solar 
system abundances of nuclei between Ba and Pb that are exclu¬ 
sively produced by the s-process. These are s-only nuclei such as 
'^'^Gd, '^'^Dy, and '^'’Yb, which typically contribute only a few 
percent to the total abundance of the eleme nt in the solar system 
(lArlandini et al.llT999l: iBisterzo et al.ll20ll1) . 

Eurther observational constraints are available for a num¬ 
ber of CEMP stars, which have metallicity roughly a factor of 
ten lower than the post-AGB stars considered here and are be¬ 
lieved to have accreted s-process material from a more mas- 
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Fig. 3. [X/Fe] ratios as a function of the atomic number Z predicted at the stellar surface by the models where one pocket 
is introduced at the end of the TDU episode after the 12* TP as compared to the observations of each of the four post-AGB 
stars indicated in the panels. We normalise both model predictio ns and observations to meteoritic rather than photospheric solar 
ab undances for the elemen ts heavier than Fe (lAsplund et al.l20()^ . Notably, this modifies the observed [Pb/Fe] upper limits reported 
by De_Smedtet^ d2014) by -0.3 dex. The black line represents pocket_case 1, which produces a similar pattern as that presented 
by 


De Smedt et al 


d2014 ) with a corresponding high Pb abundance. pocket_case2 is shown as a black dashed line. The colored lines 


represent pocket_case3 (red) and pocket_case4 (blue). (The color version of this figure is available in the online journal.) 



Fig. 4. Same as Fig. |3]but for the PIE models. The black line represents PIE_casel (Table |2]i and the two colored lines represent 
PIE_case2 (red) and PIE_case3 (blue). (The color version of this figure is available in the online journal.) 


5 
























































































M. Lugaro et al.: Post-AGB stars in the MCs and neutron-capture processes in AGB stars (RN) 


sive binary compani on that evolved through the AGB phase. 
iBisterzo et al.l (l2012h compared CEMP-i stars (showing en¬ 
hancements in the typical s-process element Ba) to AGB s- 
process models. They also included initial enhancements up to 
2 dex in the rapid neutron-capture process (the r-process) el¬ 
ements to match the Eu abundances of CEMP-i/r stars (show¬ 
ing enhancements in both Ba and the typical r-process ele¬ 
ment Eu). The s-process models give a reasonable fit to the 
observations of all CEMP-s stars for which at least one ele¬ 
ment between Eu and Pb is observed (CS 22964-161, CS 22880- 
074, CS 29513-032, CS 22942-019, CS 30301-015, CS 30322- 
023, and HD 196944). Pour out of nine CEMP-i/r stars with 
observations for the elements between Eu and Pb (CS 22898- 
027, CS29497-030, HE0338-3945, and HE2148-1247) can 
also generally be matched by the models. The remaining five 
CEMP-s/r stars with observations for several elements be¬ 
tween Dy and Pb (SDSS J1349-0229, CS 31062-050, LP625- 
44, SDSSJ0912H-0216, and HD 209621), instead, cannot be 
matched by the models because, as for the post-AGB stars dis¬ 
cussed here, the observed abundances are too high even if the 
initial r-process abundances are enhanced. 

As suggested for CEMP-s/r stars dLugaro et al.ll201^ . an 
s/r neutron-capture process intermediate between the s and the 
r process may have shaped the abundance pattern in the post- 
AGB stars. There are no detailed models of this process yet for 
low-metallicity AGB stars, though PIE events h ave been identi- 
hed as a a possible physi cal site (the i process. FCowan & Rosel 
Il977t iHerwig et al.ll201 lb . Overshoot leading to PIEs is uncer¬ 
tain, but known to be favoured in AGB stars of mass ~ 1 Mn and 


rFe/Hl<2 (Fui 

imoto et al. 19901: 

CamDbell & Lattanziol 120081 

ICristallo et alJ 

20091 Lugaro et al. 

2012h with recent 3D mod- 

els supporting these results (IStanc 

iffe et alJl201 lllHerwig et alJ 


12014 IWoodward et ahllloisll . We stress that one-dimensional 


hydrostatic models of PIE events with artificial proton distri¬ 
butions like those presented here are not expected to describe 
the ingestion and the mixing correctly and this may have a 
strong effect on the resulting abundance patterns. However, our 
ID models are still informative. Por example, we expect that 
the overabundance of Rb found in our PIE models may be a 
prominent feature a lso in more accurate, e.g., 3D, PIE models 
(iHerwig et al.ll20lTh . Purthermore, the dilution of intershell ma¬ 
terial into the envelope required to match the observations cannot 
be substantially different from that we have found here; and to 
avoid overproduction of Pb, also the amount of protons ingested 
during a PIE cannot be substantially different. Por example, tak¬ 
ing J053250 and comparing it with the model predictions from 
the PIE model that ingested 6 x 10“^ Mq of protons, we could 
match the observations by converting roughly 2/3 of the pre¬ 
dicted Pb abundance into abundances of the elements between 
Eu and Pb. In other words, a match to the observations would 
require the same amount of free neutrons that we have in our 
models, but, distributed differently among the elements between 
Eu and Pb. This may be possible if the neutron density was a few 
orders of magnitude higher than in our models and the path of 
neutron captures was shifted further away from the valley of fi 
stability, as in the i process. 

The fact that all four post-AGB stars of low mass considered 
here show similar abundance patterns in the elements heavier 
than Pe suggests that the i process may be a common occur¬ 
rence in low-mass AGB stars up to [Pe/H]~l. Since stars in this 
mass range are common this would have important implications 
for the stellar yields that drive the chemical evolution of stellar 
clusters and galaxies. Por example, these long-lived low-mass 
stars may be the /-process source required to match observa¬ 


tions of [Ba/Fe] and [B a/La] overabundances in open clusters 
dMishenina et alJ[20L5h . If the i process is confirmed to be re¬ 
sponsible for the abundances observed in CEMP-^/r and post- 
AGB stars similarities and differences in the neutron-capture pat¬ 
tern of the two groups, which sample different metallicity ranges, 
will provide fundamental constraints to pin down its metallicity 
dependence and its impact on the chemical evolution of stellar 
systems. 
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